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Abstract 

We present a comprehensive investigation of nonequilibrium effects 
and self heating in single electron transfer devices based primarily on 
the Coulomb blockade effect. During an electron trapping process, a 
hot electron may be deposited in a quantum dot or metal island, with 
an extra energy usually on the order of the Coulomb charging energy, 
which is much higher than the temperature in typical experiments. 
The hot electron may relax through three channels: tunneling back 
and forth to the feeding lead (or island), emitting phonons, and excit- 
ing background electrons. Depending on the magnitudes of the rates in 
the latter two channels relative to the device operation frequency and 
to each other, the system may be in one of three different regimes: 
equilibrium, non-equilibrium, and self heating (partial equilibrium). 
In the quilibrium regime, a hot electron fully gives up its energy to 
phonons within a pump cycle. In the nonequilibrium regime, the re- 
laxation is via tunneling with a distribution of characteristic rates; 
the approach to equilibrium goes like a power law of time (frequency) 
instead of an exponential. This channel is plagued completely in the 
continuum limit of the single electron levels. In the self heating regime, 
the hot electron thermalizes quickly with background electrons, whose 
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temperature is elevated above the lattice temperature Ti. We have 
calculated the coefficient in the well known law of energy dissipa- 
tion rate, and compared the results to experimental values for alu- 
minum and copper islands and for a two dimensional semiconductor 
quantum dot. Moreover, we have obtained different scaling relations 
between the electron temperature and operation frequency and device 
size for various types of devices. 

PACS Indices: 73.40. Gk, 06.20.Hq, 72.10.Di, 72.15.Lh, 73.20.Dx 
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1 Introduction 



Electron turnstiles and pumps have been made using metal islands [1-4] 
and semiconductor quantum dots [5,6] based on the Coulomb blockade effect. 
Other types of single electron transfer devices have also been proposed [7-9]. 
These devices have promised great potential in applications to high precision 
metrology (charge or current standard) and computer technology (memory 
and logic circuits) [l^. To date, several papers have been devoted to the 
study of possible precision of such devices, with various sources of error 
identified, such as thermal, cotunneling, nonadiabatic, and noise effects [11- 
16]. 

In this article, we present a comprehensive investigation of non-equilibrium 
effects and self heating in single electron transfer devices which are primarily 
based on the Coulomb blockade effect. In such devices, the electron tem- 
perature is very low compared to the Coulomb charging energy, but may 
exceed the level spacing in the excitation spectrum. Devices based on metal 
islands clearly lie in this regime, and this is also true for currently available 
quantum dot turnstiles or pumps, where excitation level spacing is typically 
several times smaller than the Coulomb charging energy. For ease of refer- 
ence, we will use the term Coulomb island for a metal island or quantum 
dot in such a regime. Some highlights of this investigation have previously 
been published [jl3|, ITSf. The primary purpose of this article is to provide 



a complete picture of the different regimes, with comprehensive arguments 
and calculations included. A number of new results will also be presented. 

We first examine a quantum dot electron turnstile to identify a major 
source of non-equilibrium and self heating: the deposition of hot electrons 
in a Coulomb island |]18|. This process is common to all Coulomb blockade 



devices. Figure 1 shows the process of an electron tunneling from a lead into 
the quantum dot, which occurs when the barrier between them is lowered 
[|, §]. The bias is large enough to overcome the coulomb charging energy 
(U) for the addition of one electron into the dot, but not so large as to allow 
a second electron to enter the dot. To minimize the thermal probability of 
adding no electron or more than one electron into the dot, the Fermi energy 
of the feeding lead is typically placed at about U/2 higher than the lowest 
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available state in the dot before tunneling. There are then several available 
states into which an electron can tunnel from the lead. It is imaginable that 
the electron may actually tunnel into a state of about U/2 higher in energy 
than the lowest available state in the dot, the probability of which is further 
enhanced by the fact that the tunneling rate increases rapidly with energy. 

There is a similar non-equilibrium process during the draining part of a 
cycle. The Fermi energy of the draining lead is usually biased many levels 
below the highest occupied one. When an electron tunnels out to the draining 
lead, it may come from a level below the highest occupied one. As a result, 
a hot hole may be stuck in the dot or the metal island for a long time due to 
the weak coupling to the lead and to the phonons and background electrons 
at low temperatures. Since the hot hole effect is entirely similar to the hot 
electron effect, we will not consider it explicitly for most of this article. 

We then consider an electron turnstile based on small metal islands as 
shown in Figure 2. Unlike the quantum dot device, the tunneling barriers 
cannot be varied, and electron transfer is purely controlled by changing the 
various potential biases in the circuit. To trap an electron in the middle 
island from the left lead, a bias larger than e/C is applied between them, 
where C is the capacitance of the tunnel junctions. One possibility is that 
an electron will enter from the lead into the left island, causing the Fermi 
energy of the left island to be ?7 = e^/2C higher than the lowest available 
state in the middle island. As a result, a hot electron will be deposited into 
the middle island from the left island. An alternative possibility is that an 
electron tunnels first from the left island to the middle island followed by the 
deposition of a hot electron in the left island from the left lead. Finally, there 
is the possibility of co-tunneling, which may also result in a hot electron in 
the middle island |T2| . 

The situation in an electron pump is slightly different (Figure 3). The 
tunneling through a junction occurs while the bias across it is ramped up 
through the threshold of Coulomb blockade. In the adiabatic limit, tunneling 
will occur at the threshold, with one electron deposited in a lowest available 
level in the final island. However, a hot electron will be in the final states, if 
the ramping rate is finite. On average, the extra energy of the hot electron 
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increases as a square root of the ramping rate 



Therefore, a non-equihbrium situation can easily arise in a Coulomb island, 
unless the hot electron can relax to the lower available states in a time short 
compared to the operation cycle of the device. There are three ways of hot 
electron relaxation: (1) tunneling back and forth, (2) emitting phonons, and 
(3) thermalizing with other electrons. In the first process, electrons tunnel 
back and forth between the Coulomb island and the feeding lead (or island), 
eventually bringing the two regions into mutual equilibrium. However, such 
a relaxation process can be extremely slow at low temperatures, because 
the rate of tunneling in a second electron is severely reduced by the Coulomb 
blockade, and that of tunneling out the hot electron is very effectively blocked 
by the lack of holes in the Fermi sea of the feeding lead (island). In the sec- 
ond process, the hot electron simply emits a phonon and makes a transition 
to a lower state. At typical energies of order 0.3 meV or lower, only acoustic 
phonons are emitted, and the rate of emission is significantly limited by the 
weakness of electron-phonon coupling and the lack of phase space. Finally, 
in the third process, the hot electron relaxes at the expense of exciting other 
electrons already present in the island. The hot electron may reach equilib- 
rium with others in the island by electron-electron scattering, but the island 
as a whole is heated. 

There are three regimes of operation according to how fast the rates of 
the relaxation processes are compared to one another and to the operation 
frequency of a device. These are illustrated in Figure 4 as different regions in 
the jee-leph parameter space, where 'jee is the characteristic rate of the hot 
electron relaxation due to interaction with other electrons, and 'jeph is the 
rate due to emitting phonons. When both rates are slow compared to the 
operation frequency u of the device, we have the non-equilibrium situation. 
When the phonon emitting rate 7eph is larger than both the operation fre- 
quency u of the device and the thermalizing rate 7ee, we have the equilibrium 
situation. On the other hand, when the thermalizing rate jee is larger than 
both z/ and 'jeph, we have the self heating (partial equilibrium) regime. 

The article is organized according to the different regimes. The non- 
equilibrium regime is studied in detail in Section 2, where only the tunneling 
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process is taken into account. One will see that the relaxation of a hot 
electron is characterized by a set of distributed rates, the lower spectrum 
of which are severely frozen at low temperatures. As a new result, we will 
establish a power law frequency dependence of energy relaxation of the hot 
electron when the level spacing in the Coulomb island becomes smaller than 
the temperature. Some details of this section are included in the Appendix. 

In order to delimit the equilibrium regime, we calculate in Section 3 phonon 
induced interlevel transition rates, and study how the relaxation rates found 
in the previous section arc affected by the phonons under experimental con- 
ditions. Discrete level effects are emphasized for applications to quantum 
dots, and the continuum Fermi liquid model is considered for metal islands. 

Section 4 is devoted to the study of the self heating regime. In the first part 
of this section, is estimated for typical quantum dots and metal islands. In 
the second and third parts, the rates of energy dissipation from thermalized 
electrons to phonons are obtained for a metal island and a quantum dot, 
respectively. The coefficient in the law of energy dissipation is calculated 
for aluminum and copper islands and for a two dimensional quantum dot, 
and these new results are compared to experimental values. In the last part 
of Section 4, we show how the electron temperature and thermal error scale 
with the frequency and size for different devices. Finally, we summarize in 
Section 5. 

2 Hot Electron Relaxation Via Tunneling 

In this section, we simulate numerically the action of an electron turnstile, 
and show analytically how a hot electron can appear and relax in a Coulomb 
island during the trapping part of an operating cycle. As mentioned in the 
introduction, we will concentrate our attention on the tunneling processes, 
with studies of the effects of phonons and background electron excitations 
deferred to the next two sections. Our discussions will be mostly focused on 
a quantum dot device with discrete levels, but we will also consider the case 
where the level spacing is small compared with temperature in the leads, 
and compare the results with existing theory about tunneling between metal 
islands. 
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2.1 Modeling the Quantum Dot Electron Turnstile 



Our model consists of a quantum dot separated from two leads by tunable 
energy barriers. The leads are represented by two electron seas with Fermi 
energies Ej and respectively. We will consider a set of available single 
electron levels with spacing A in the dot. We place Ej in the middle of the 
gap above the [yjth level (counted from below), and Ej^ at U/2 below the 
lowest available level, where U is the Coulomb charging energy in the dot 
[26, 20 1 . Corresponding to current experimental conditions, we will assume 



hu -C ksT <^ U and hT << A, where u is the operating frequency and fiT 
the level width due to tunneling. Then semiclassical rate equations can be 
used to describe the tunneling processes: 



dt 



(1) 



l=L,R 



where = [l + exp( ° j-^ )]'^ is the Fermi-Dirac distribution function in the 



leads. Pa is the level occupation probability in the quantum dot, and is the 
coupling of the dot to the leads. The factor (1 — Y^a' Pa') in the tunnehng-in 
term reflects the Coulomb blockade effect in the large 



u 



2kBT 



limit. The rate 



equations have been applied to many different situations, such as the study 
of conductance oscillations in a quantum dot device [1^, 0, single electron 
tunneling through metal junctions |2l|, charge transfer scattering of atoms 
or ions from a metal surface ^ and non-equilibrium effects in a turnstile 
dot device |13]. In the Appendix, we give a brief summary of a study on the 



non-Markovian correction to the above rate equations E^, E3, showing that 



the semiclassical rate equations are accurate enough even for the purpose of 
examining exceedingly small errors of electron trapping and transfer. 



The couplings between the states in the quantum dot and the leads are 
modeled as follows. To simulate the sensitive energy dependence of the cou- 
plings, we take g^it) = g^j{t)a"~"^ , where a > 1 is a constant, and «/ is 
the label of the state immediately below Ej ||2^. To move one electron 



from lead L to lead R per cycle, we adopt the time-periodic couplings as 
g^{t) = rj^e-'^(^+™'^('^*)) and gf{t) = r{?e-^(^-™"(^*)), where S characterizes 
the variation amplitude of the coupling strengths. It can be seen that in the 
first half of the cycle, the coupling to the right is negligible and only those 
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electrons in the left lead is allowed to move into the dot; while in the second 
half, the left lead is nearly blocked and the trapped electron moves to the 
right via the large coupling to the right lead. For simplicity, we assume that 
the energy levels do not change with time. 

2.2 Trapping Probability and Error 

The basic operation principles of all single electron transfer devices are 
based on an adiabatic analysis appropriate for very low operation frequencies. 
In the adiabatic limit, the time dependence of the occupation probability of 
the dot, -P(t), is found by letting the right hand side of Eq.(l) vanish, yielding 
P{t) = S[l + Sr' with S = Ec. [(1 - fa) + gf{t)/gj{t)r\ where we have 

dropped the terms of order under the assumption that e ^''b'^ is as small 
as desired. At t = 0, where gf/gf takes the maximum value of order e^^, 
S and therefore P(t) take a minimum value of order e~^^. This is due to 
leakage through the left barrier which is high but finite at t=0. Since we are 
only interested in non-equilibrium effects due to finite operation frequencies, 
we assume that the leakage error is also as small as desired. On the other 
hand, at t = tt/cj, where gf/gf takes the minimum value of e~'^^ << 1, 5* is 
of order of (e'^/^'^s^ + e~'^^)~^, and P(t) approaches unity with a deficiency 
of the order of (e"'^/^'^^-^ + e"^''), which is to be ignored according to the 
above discussion. 

We have numerically solved for the occupation probability P{t) and tun- 
neling currents as functions of time for finite operating frequencies. Shown 
in Figure 5 is the case for u = 20 MHz, with realistic parameters taken from 
Ref.|, I as (5 = 10, ^ = 4, a = 1.5, = = 10^° /s and = 10. 
The time interval shown is taken after the first cycle so that the memory 
effect of the initial condition is absent. The general trend in P{t) follows 
the adiabatic result, but the minimum point to and the maximum point tm 
are shifted by 0.27r above out = and by 0.327r above ut = vr, respectively. 
The tunneling currents from the left lead and to the right lead are given by 
Ii = ±(7i(t)[(l-Ea' Pa')fL-Pail-fL)l where l = L,R. They are shown in 
Figure 5 to aid our understanding of the electron trapping and expelling pro- 
cesses. It is noticed that when P(t) rises or drops down to | at ti/2 = 0.72% 
or 1.727r, the current II or J/j reaches its maximum. In the interval of (to, tm) 
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the unwanted leakage is negligibly small (< 8 x 10~^uj), this corresponds 
to the trapping process; while in the interval {tm, h + l/i^) the leakage II is 
neghgibly small (< 9.8 x 10~^a;), this corresponds to the expelling process. 
The probability that an electron is transferred from the left to the right lead 
in a time cycle is Jta^^^'^ ~ ^t^ ~ Max{P) — Min{P), where in the 
first approximation we have neglected II in the second half of the cycle, and 
in the second approximation we have neglected in the first half. Also, 
we found that Min{P) ~ 10~^. Therefore, up to terms of order 10~^, the 
electron transfer probability is the same as the electron trapping probability 
Max{P) in the first half cycle. 

We would like to know how the trapping error varies as a function of the 
operation frequency. In Figure 6, the trapping error (1 — Max{P)) (open 
diamonds) is plotted against 1/v on log- log scales. Except for the very lower 
end, the curve also represents the transfer error according to the above dis- 
cussions. The general decreasing trend of the curve appears like a power law 
decay, but the downward curving segments are actually closely represented 
by exponentials. We explain these points in the following analysis. 

2.3 Distributed Relaxation Rates 

To explain the frequency dependence of the trapping error, we closely 
examine the time interval | < i < (sec Figure 5), in which the coupling 
to the right lead can be ignored. We can then write 

a' 

We first find the equilibrium (adiabatic) solution P^^ of the equation by let- 

ting the right hand side equal zero, yielding: P^^ = e ''b'^ (i -|- S)~^, where 

S — J2a 6 '"^^ ■ The equilibrium occupation numbers obey the Boltzmann's 
distribution, and the sum J2a approaches unity with a negligible correc- 

-u 

tion ~ e'^'^BT _ 

Next, we consider the actual probability distribution, and see how it differs 
from equilibrium. For the time interval | < t < ti/2, we may assume that 
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the occupation of each level is negligible, so that the term proportional to 
Pa on the right hand side of Eq.(2) may be dropped. We then find that 
Pa{t) is proportional to a°'~°'f f^, with the a-independent coefficient equal 
to /| dt'g^{t')[l — P{t')\, where we have neglected the initial value Pa{\)- 
The coefficient can be fixed to the more explicit form, P{t)/ 0"""^/^ > by 
using the condition that P^ should sum up to P. It is clearly seen that the 
actual probability distribution in this time interval is far from the equilibrium 
Boltzmann distribution. In fact, it is peaked at the level just below the Fermi 
energy, corresponding to the tendency of depositing a hot electron there. This 
distribution is well reproduced in the numerical calculation even at i = ^1/2- 

To see how this distribution relaxes to equilibrium in the interval, ti/2 < 
t < tm, have performed an eigenvalue analysis. The deviations, APa{t), 
from the equilibrium occupation numbers may be written as 

AP„(t)=E^ie-*'^'U (3) 

3 

where the A's and A's are eigenvectors and eigenvalues (labeled by j) of the 
following equation 

XAa = a^-'^mHAa')!^ - Aa{l - f^)]. (4) 

a' 

The eigenvectors are only determined up to normalization factors, but they 
can be fixed by using the values of APa{ti/2), which are known from the 

-u 

discussions in the preceding two paragraphs. Up to terms of order e^'^s'^ , the 
trapping error l—Max{P) is just the negative of J2a ^Pa{tm), and is given by 
EjB^e-"^/", where = - Ea^i and Uj = \jq{t^)u ^ \jV^f/V^, which 
is frequency-independent. The trapping error calculated this way is plotted 
in Figure 6 as the solid curve, which fits the result of two-lead numerical 
simulation, represented by the open diamonds, very closely for the entire 
range of several decades of data. 

The eigenvalues are listed in Table 1 for the parameters N = 10, a = 1.5, 
and A/ZcfiT = 4, which are used for Figure 6. Note the huge differences in 
the magnitude scales of the eigenvalues. At 5 = 10 and — 10^° /s as used 
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in the figure, they give rise to the following series of characteristic relaxation 
rates: z/i = 0.13 x 10^ /s, i/g = 14 x 10^ /s, 1/3 = 1.2 x 10^ /s, z/4 = 96 x 10^ 
/s, and etc. The three exponential segments seen in the figure are attributed 
to V2 to 1/4. One can also predict that at very low frequencies, the exponent 
should be given by v^. The power- law- like overall trend of the curve is a 
result of distributed exponentials. 

We now present an analysis of the eigenvalues in order to see their physical 
origin. If we solve for A^, from the eigenvalue equation, and make a sum- 
mation over a, we find that A satisfies the equation Z^jLix^ ~ where 
dj — a^~"-f{\ — fj") and bj = a^'^'f f^. It can be seen that \j lies in be- 
tween dj and dj^i for j = 1, ...A^ — 1, and that A^v lies above d^^. All the 
eigenvalues are positive as required. Because of the Pauli exclusion factor 
(1 — /j^), the dj's with j < af are exponentially small at low temperatures, 
which, together with the fact that the tunneling rates increases rapidly with 
level energy (a > 1), explains why the lower eigenvalues become so small. 
Quantitatively, each of the roots can be approximately found by keeping only 

the neighboring terms in the summation series, yielding \j = ^ i+a~"^^^^^ 

for 1 < j < a/, and Xj ~ a-' "^ for those j > af. These rough estimates 
agree well with the numerical calculation. 

We have also done a similar calculation for = 2, and found that the 
characteristic frequency scales arc raised to ui = 0.15 x 10^ /s, z/2 = 2 x 10^ 
/s, z/3 = 24 X 10^ /s, z/4 = 227 x 10^ /s, and etc.. The qualitative behavior 
may be understood by the fact that the Pauli exclusion effect is less strong 
at a higher temperature. Finally, we note that the characteristic frequency 
scales may also be increased (in our favor) by increasing Fj and lowering 5, 
but one must pay the price of increasing the unwanted leakage. 

2.4 Power Law of Hot Electron Relaxation 

We mentioned in subsection 2.2 that the overall trend of the trapping 
error looks like a power law of frequency. To explain this, we consider in 
this subsection the case of small level spacings compared to the temperature. 
We expect that the exponential segments in Figure 6 will become invisible, 
and that the frequency dependence of the trapping error will become a pure 
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power law. Such a study is particularly relevant to devices based on metal 
islands [1-4], where the single particle levels form a continuum. 



The major tool of our analysis is the method of Laplace transformation. 
We first rewrite Eq.(2) as 



dR 



a 



a"-"/[(l-^P.O/a-^a(l-/.^)], (5) 
a' 

where s = /^/(2(^) g^{t')dt' . Performing a Laplace transformation with respect 
to this new 'time' variable, we obtain that AP„ = a"^"^[(l/A — P)/^ — Pa(l~ 
fa)]-i where A is the Laplace variable and ~ denotes transformed quantities. 
We have used the initial condition that the probabihties are zero. Solving 
for Pa in terms of P and summing over a yields an algebraic equation for 

P, which has the solution P = ;xT+5' where S — J2a x+(i'-"fL)a°'~°'f ' ^^^^ 
then obtained from the standard formula of inverse Laplace transformation: 
^e^^P(A), where the contour surrounds the singularities of P(A). 

The pole at A = yields the equilibrium probability Peg = x+5(o) i "which 
is the same as before. The other singularities are poles lying on the negative 
part of the real axis, which become a branch cut along the segment (—1, —ho) 
in the limit of densely distributed levels, where ho — 1 — fo « 1 correspond- 
ing to the lowest available level. In the asymptotic region 1 << s « 1/ho, 
we may evaluate the trapping error analytically as 

P - Pis) ^ ^M^^^lO^l^s-^ (6) 

where 77 = (1 + lna)~^, A is the level spacing, and r() is the standard F 
function. The case of a = 1 (r^ = 1) needs special consideration, where the 
second factor on the right hand side of the above equation should be replaced 
by the logarithm (Ins/io)^^; however, in general, a > 1 and rj < 1 holds since 
the coupling to a dot or a metal island decreases as the level energy drops. 



Since at t — tm we have s — / {i'\/2-k5), the above asymptotic range 
corresponds to voho « v « i^o, where vo = Tj/\/27r5. The above result 
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says that in this range the trapping error goes as a sub-hnear power law 
of frequency v. Corresponding to the parameters used in Figure 6 we have 
7] = 0.9 and uq = 1.26 GHz. The best fit to Figure 6 yields a power of 
T] = 0.8, which is quite close to the above prediction, considering the fact 
that the parameter A/{kBT) = 4 used for the figure is far from the limit of 
dense levels. 

For an electron turnstile using metal islands, the tunneling rates do not 
vary with time {6 = 0), and we should have s ~ r^/{2i') at t = t^, yielding 
the same power law behavior of the trapping error in the asymptotic range 
hoTj << V << rj. However, the electron level density in a metal island is 
much larger than that in a quantum dot, and correspondingly the tunneling 
rate Tj for a single level may be very small. It is therefore more appropriate 
to consider the regime of u » Tj, where we may drop the second term on 
the right hand side of Eq.(2). Then a single simple equation for P is obtained 
by summing over a, with the result 

r/P 

— = (i-p)rjx:«"-'^^/^ (7) 

Therefore, the trapping error relaxes exponentially with a single rate 7 given 
by the quantity multiplying the factor (1 — P) on the right hand side of 
the above equation. The trapping error at t = tm goes exponentially as 
^-iT"//u^ For a = 1 and at low temperatures, this rate is the well known result 
7 = [Ej' — EQ)/{e^RT)i where is the Fermi energy of the feeding lead 
(island), Eq is the energy of the lowest available level in the trapping island, 
and Rt = A/(rje^) is the tunneling resistance. 

However, the exponential relaxation of the trapping error does not mean 
that the distribution of probabilities on the levels also relax exponentially 
to equilibrium. In fact, with the known time dependence of P substituted 
back into Eq.(|^, we found that Pa has the long time limit of Tja°'~°'f f^/j, 
which is far from the Boltzmann distribution at equilibrium with the feeding 
lead. Indeed, the extra energy J2aPa{Ea — Eq) of the trapped particle is 
found to be on the order of Ej — Eq instead of fc^T at equilibrium. This 
failure of reaching equilibrium in the continuum limit can be understood. 
The second term on the right hand side of Eq.(2) describes the process of 
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tunneling out of the trapped particle, a necessary step for the redistribution 
of the probabilities on the levels, but we have dropped it due to the vanishing 
of in the continuum limit. In reality, such a term does exist although it 
is very small, and equilibrium will be reached after a very long time. We 
will not look further into this, because this time scale can be much longer 
than an operation cycle of the device, or the time scales of inelastic processes 
involving phonons or background electrons. 



3 Hot Electron Relaxation Via Phonon Emis- 
sion 

The above discussions have not taken into account interlevel transitions, 
which may facilitate the trapped electrons to attain equilibrium within the 
dot as well as with the ambient. In this section, we will estimate the transition 
rates due to phonon emission. Both deformation potential and piezoelectric 
coupling will be considered. Particular attention will be given to the case 
of a quantum dot, in which the discreteness of the energy levels plays an 



important role ^7j. For completeness, we have also considered the contin- 
uum model, which is appropriate for metal islands and for sufficiently large 
quantum dots. 



3.1 Deformation Potential Coupling 

Following common practice ||2^, we write the deformation potential cou- 
pling as: Hd = J^jk^jk^jk, where Ej^ is the deformation potential tensor, 
and ejk is the symmetric strain tensor ^ + |^] [20] . The lattice displace- 
ment of the substrate may be expanded in terms of the normal modes as 
u = Sg,A[ 2A/c;|g| ]^ + (^-q,\)]: where M is the total mass of the 
substrate, q is the phonon speed, and agx and the annihilation and 

creation operators of a phonon with wave vector q and in mode A, respec- 
tively. We have treated the lattice points Rj by a continuous variable R, 
which is appropriate in the long wavelength limit. The polarization vector 

is required to be real and satisfy = — ^-g. 
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For an electron initially in state \i > with energy Ei and finally in state 
1/ > with energy Ei+SE, the transition rate due to the deformation potential 
is obtained using the Fermi's golden rule as: 

' Mci^ q 



SE 



where hb = {ub + | ± 2) with ub^SE) = [e^s^ — 1]~ , the positive sign is 
for 5E < and negative for SE > 0. The stress tensor Ejk may be reduced 
to a diagonal form using its principal axes, and the sum involving Ejk is 
simplified to | Z^ja -^ji(^q,A'?i)P- ^ spherically symmetric conduction band 
valley, Ejj = D. Then, Eq.(6) becomes 

rS? = ^^^E^^l < ¥'-V > \'5{hcM\-5E)nB. (9) 

For a nonspherical band valley, the stress tensor components may be different 
from each other, but we may still use Eq.(||) by replacing the sum in Eq.(||) 
by its average over the direction of q. We see from Eq.(^ that only the 
longitudinal phonons, for which ^ ■ = Q'^? are involved in the transition. 
Replacing the summation J2q by the integral x / d^q, where V is the 
volume of the substrate, we derive that 



where p is the mass density of the substrate, Mif = \ < i\e^'^'^\f > p, and Mjj 
is the average of transition matrix element M^f over all possible directions 
of q for a given energy difference 6E. For GaAs, the deformation potential 
D = 8.6 eV, phonon speed q = 5136 ms~\ and mass density p = 5.3079 
g/cm^ [^, then Eq.(lO) becomes 

rS? = 5.3xlOV^x(^^fMn.. (11) 

3.2 Piezoelectric Coupling 

The piezoelectric potential (pp satisfies the equation 
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where ptjk are the piezoelectric constants. For detailed discussion of the 



piezoelectric properties of GaAs crystal, the reader may consult Ref 
Since ejk can be expanded in terms of the normal modes as 

the piezoelectric potential energy is obtained from Eq. (9) as 

= eE[lTir^]M2.Ae^^^-%A + C.C.], (13) 

where Eqx = Eijk ^^ixPijk- 

Applying the Fermi's golden rule again, we obtain the transition rates due 
to the piezoelectric coupling as: 

r!? = < eo6.E^ > I < ^le^^^-^V > mhcM\-5E). (14) 

eoerV q q x P'^i 

Replacing < tQtr J2x > by the piezoelectric constant p'^ and the 

pCj 

summation over q by its corresponding integral, we obtain 

A/? 

= 2.7xl0iVix— -M.ffiB. (16) 

O.lmeV 

In the second equality, we have used the values of the piezoelectric constant 
p = 0.052 and the relative dielectric constant Er = 12.91 for GaAs p9|. 



The above result does not contain the effect of screening from the electrons 
inside and outside the quantum dot, and may be used as a first approxima- 
tion when the quantum dot contains only a few electrons. However, for a 
quantum dot of few thousand A, hundreds of electrons are present |^, |], 
which considerably reduce the piezoelectric potential. In such a case we may 
regard the electrons in the quantum dot as a 2DEG, for which the screening 
factor for the transition rate is given by (asgy) , where 
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is the screening length of a 2DEG in a GaAs heterostructure, and q\\ is the 
phonon wave number in the plane of the 2DEG. For a rough estimate, we 
may take q\\ as mt/ where n is the number of levels between the initial and 
final states. Taking a level spacing of 0.03 meV, and dot size L = 4300 A, 
the screening factor may be written as 0.014(5i?/0.1 meV)^ ||. 



3.3 The Transition Matrix Element 

We now consider the transition matrix element M. Some detailed discus- 



sion of M for different device geometry in a 2DEG has been given in Ref. 
To understand the q dependence of M, we consider, for instance, the case 
of a square dot with the transition (fc^., ky, k^) = {IOtc/L^, lOn/Ly, vr/L^)) —>■ 
(IOtt/L^, 97r/Ly, tt/L^), where the fc's are the wave numbers of the electron 
state sm{kxx) sin{kyy) sin{kzz). Using Eq.(4) of Bockelmann and Bastard 
we obtain, after some mathematical manipulations, the transition ma- 
trix element as M = '^'"^ ^nl-%^2^2 ? where we have introduced the re- 
duced phonon wave numbers: = q^Lx/l, Qy = qyLy/2, and = qzLz/2. 



Now, it is instructive to give an estimate of the phonon wave length, using 
the law of energy conservation: A = hci/6E. For 0.03 meV < 6E < 1 meV, 
we find A ~230-7000 A. Since the thickness of a 2DEG is only ~ 100 A, 
we can set M^(g^) = 1. Then, the average of the transition matrix element 
over the directions of the phonon wave vector are found as: M~|(^)'^Q^ for 
Q<Cl and ~ for Q S> 1, where Q = qL/2, assuming = Ly = L. 
For transition energies 6E > 0.03 meV (level spacing) and a dot size of 
L = 4300 A, we find Q > 2. Employing the formula for large Q, we obtain 
M = 2xl0-2(^^)-2. This gives 

for deformation potential coupling, and 



Tt} = 7.6 X 10« X 5^«B (18) 

for piezoelectric coupling (with screening). It is seen that both deformation 
potential and piezoelectric couplings give a transition rate of similar form, 
but the latter is more effective for GaAs. 
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Some comments are in order. First, although we have considered a par- 
ticular pair of initial and final states, the results also approximately hold for 
other transitions in the dot. Second, the large Q form of the matrix element 
can also be obtained between a pair of plane wave states. Comparison with a 
continuum model will be discussed in subsection 3.5. Third, when the wave 
length of the emitted phonons are longer than the size of the dot, the matrix 
element in the small Q limit should be used, for which the transition rate 
Fj/ ~ {^EY. The scaling of Q with the dot size is somewhat tricky. For 
fixed 5E^ Q decreases linearly with the size. However, since the level spacing 
scales as , Q increases as with decreasing the size, if the number of 
levels between the initial and final ones is fixed. 



3.4 Effects on the Distributed Relaxation Rates 

To study the effect of phonon induced intradot transitions, we add the 
rates determined by Eq. ([T8|) to the right hand side of Eq.(2): 

dP ^ 

"-^ a' j=l 

where Tj^ = F-j''. Since the phonon induced transitions leave the Boltzmann 
distribution invariant, the equilibrium solution remains the same as before. 
However, the non-equilibrium behavior is changed substantially. To get an 
idea of how the intradot transition rates modify the trapping process, we 
approximate g^if) by Fj for \ut — ti\< and by otherwise. 

The trapping error has the form of Eq.(3) but with different eigenvalues 
{Aj} and eigenvectors. As we can see from Table 1, dramatic changes occur 
in the slower modes represented by the first four columns. In particular, the 
lowest characteristic rate is now raised to vi = 10~^rj/V27T6 ^ 12 x 10^ 
/s. This is about the same as the phonon-induced rate for an electron to 
relax from the level just below the Fermi energy of the left lead to the lowest 
available level in the dot, which is Fj/ at 6E = 5 x 0.03 meV. We have also 
solved a similar case with = 2, and find that the lowest characteristic 
rate is raised to z/i = 8.4 x 10^ /s. Thus, electron phonon interaction helps 
to thermalize hot electrons in a Coulomb island. 
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3.5 The Continuum Limit 



In this subsection we calculate the phonon induced decay rate of a state in the 
Fermi liquid model. This is appropriate for a metal island in which the levels 
are densely distributed. The result for a quantum dot in the continuum limit 
is also useful as a reference for our understanding of the results calculated 
for the discrete case. 

— * 

Consider first a metal island, in which a state k at energy 5E above the 
Fermi sea is to decay into other states via phonon emission. The rate at zero 
temperature is given by the Fermi golden rule as 

leph = y E W^iJ^^i -Ek + Ek,), (20) 

k' 

where q = k' — k is the phonon wave vector and ujg = qci is the phonon 
frequency. The coupling function is taken as \gq\'^ = where n is the 

electron density and V is the volume |3l|] . Replacing the sum by a k' integral 



and carrying out the radial integration, we obtain 

where ds' = 2'Kk'lsin6 d6 is a Fermi surface element, and 6E << Ef is 



bin c uc lb a reiiiii buiiace eieiiieiiL, anu uni nif 
assumed. Employing the relations kj sin 9 d9 = qdq and n = kj/{37i'^), we 
carry out the resulting q integral in the allowed range < g < 6E/{hci), 
yielding 

where m is the electron mass. The cubic power dependence on energy can 
be understood by a simple power counting: 2 from the phase space of the 
phonons, and 1 from the coupling constant. 

We now make some estimates. For aluminum, Ej = 11.7 eV and q = 6420 
m/s, we have 'jeph = 0.072 x 10^ /s at 6E = 0.1 meV. For copper, Ef = 7 eV 
and Q = 5010 m/s, we have 'jeph = 0.2 x 10^ /s at 6E = 0.1 meV. 
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Next, consider a quantum dot in the continuum limit. Using a similar 
method, we obtain the decay rate of a state due to emission of 3 dimensional 
phonons as: 

where is the phonon wave number in the direction perpendicular to the 
2DEG. Replacing the Fermi line element dl' by dq\\, and carrying out the 
angle integration in the {qz,(l\\) plane, we have 

For deformation potential coupling, we have V^ls'gP — ' yi^^^ii^S 

^ zAnpcfn Vf 

For a 2DEG in a GaAs/AlGaAs heterostructure, we may use Ej = 7 meV 
and effective mass of an electron m=0.067me, yielding Vf = 1.9 x 10^ m/s. 
Using the known parameters, we find ^eph — 1-2 x 10^ /s at 5E — 0.1 meV. 

For piezoelectric coupling, we have ^l^fg^ = Tie^ Cip^ / {2eQerq) . This is to 
be screened by a factor of {q\\asY^ which becomes \{qas)'^ after directional 
average in the (g^, q\\) plane. The decay rate then becomes 

leph = — TT^ , (26) 

which yields 'jeph = 0.6 x 10^ /s at 6E = 0.1 meV. 



Some comments are in order. First, the cubic power dependence on en- 
ergy remains true for the 2DEG, because we still used 3 dimensional phonons. 
Second, in subsection 3.3 the transition rate from a given initial state to a 
given final state was found to be ~ 5E. This is consistent with the (SE)^ law 
of decay rate found in the present subsection, because here we have summed 
over the final states which occupy a phase space of order {SE^. Third, for 
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the example considered in the previous subsections, there are only few avail- 
able final states, so the transition rate found there also corresponds roughly 
to the decay rate of the initial state. This rate is, however, considerablly 
larger than the decay rate predicted by the continuum model, illustrating 
the importance of discreteness of the energy levels. Finally, formulas derived 
in this subsection may be used for a crude estimate of the cooling rate of an 
electron gas at a temperature Te higher than the lattice temperature. Take 
an electron at energy SE ksT^ above the Fermi energy, it will decay into 
lower states with a rate ~ (kBTg)^ and with an energy transfer of order ksTe. 
The number of such electrons per unit volume (area) is of order ksTe times 
the density of states at the Fermi energy. The cooling rate is therefore pro- 
portional to and to the volume (area) of the Coulomb island. Quantitative 
calculations of the cooling rates will be done in the next section. 

4 Self Heating in a Coulomb Island 

We now bring in an additional channel of relaxation of the tunneled-in 
electron: the inelastic scattering by other electrons in a Coulomb island. Such 
scattering has two effects: (1) it limits the characteristic frequencies of the 
device to be above the typical rate 7ee of this channel, and (2) it thermalizes 
the tunneled-in electron with others, turning the extra energy into electronic 
heat. If the thermalization process is fast enough, we may assign an electron 
temperature Tg for the electron gas in the island, although it may be higher 
than the temperature of the lattice phonons in the substrate. Phonons will 
be emitted from the electron gas to reduce the temperature difference. A 
steady state is reached when this cooling rate is balanced by the heating rate 
due to tunneling- in hot electrons (and holes). 

In the following subsections, we first give a rough estimate of the hot 
electron relaxation rates in both a metal island and a quantum dot due to 
interactions with the background electrons there. We then calculate the 
coohng rates of thermahzed electron gases due to emission of phonons; the 
results for metal islands of aluminum and copper and for a quantum dot 
will be compared to experimental measurements. Finally, we show how the 
electron temperature (and the corresponding thermal error) in a steady state 
scales with the size and operation frequency for different devices. 
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4.1 Relaxation Due to Background Electrons 



For a sufficiently clean and large metal island, we may use the Fermi liquid 
formula 7ee = ^'j^p to estimate the inelastic relaxation rate of an electron 
of energy SE above the Fermi energy Ef of the metal island, where we have 
assumed that the temperature of the background electrons is much smaller 
than 6E [^. The relaxation rate is found from this as 0.2 x 10^ /s for an 



Al island with 6E =0.1 meV. 



The correction due to disorder enters in a somewhat complicated manner 

depending on the length scale Lse = \J'j^i where / is the elastic mean free 
path |3^. In order to have an idea of how big this length scale can be, we 
take 5E = 0.1 meV and / = 100 A, then Lse ~ 4000 A for aluminum or 
copper. If Lse is smaller than the linear sizes of the island in all directions, 
then an additive correction to the relaxation rate is given approximately 
by {hNfL^^)~^ ~ {SEY^'^, where Nf is the density of states at the Fermi 
energy. This correction becomes larger than the Fermi liquid result when SE 
is smaller than the level spacing in a mean free volume l^. For 6E = 0.1 
meV, the crossover occurs when I is about 80 A for aluminum or copper. 



On the other hand, if Lse is larger than the linear sizes of the island in all 
directions, then in the above formula should be replaced by the volume V 
of the island. In this regime, the correction to the relaxation rate is a constant 
in 6E and becomes larger than the Fermi liquid result at 6E = 0.1 meV when 
the volume V becomes smaller than 0.05 /im^. Such an island volume is in 
fact quite typical in Coulomb blockackade experiments. Finally, the system 
can be in neither of these regimes if Lse lis in between different length scales 
in different directions of the island. For example, for a metal island of sheet 
geometry, the correction to the relaxation rate goes as {hNfLsEd)^^ ~ 
if Lse becomes larger than the thickness d but still smaller than the other 



length scales of the island. For more details, the reader is refered to |34 



For a semiconductor quantum dot with a degenerate 2DEG, the Fermi 
liquid formula for the relaxation rate of an electron of energy SE << Ef 
reads: 



(5^)2 1 , ^ IK 



, +ln2 + -ln^l. (27) 

hEf '4: 2 5E^ ^ ^ 
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We have obtained this result by following the method of Hodge et al 
and using the two dimensionally screened Coulomb interaction. The factor 
including the logarithms is due to the two-dimensionality of the problem, and 
is numerically about 2 to 3 for 5E = 1 ~ 0.1 meV, where we have used Ef = 7 
meV for a GaAs/AlGaAs structure. The relaxation rate is then estimated as 
7ee = 3.5 X 10^ /s for 6E = 0.1 meV, which is three orders magnitude higher 
than that for a metal island due to the much smaller Fermi energy for the 
2DEG. 

Corrections due to disorder in the strictly 2D case is given by g^^^, where 



m is the effective mass [^]. (It is interesting to note that this can also 
be written in the suggestive form {hNfL'j^)~^ up to a factor of 4, where 
Nf = is the density of states in two dimensions.) For 5E = 0.1 meV 
and Ef = 7 meV, we find that the correction becomes larger than the Fermi 
liquid result when / < 1700 A, where we have taken m to be 0.067 times the 
electron mass. 

4.2 Cooling Rate of a Metal Island 

We assume that the electron gas in a metal island has a well defined tem- 
perature Tg, which may be different from the lattice temperature T;. When 
a phonon of momentum q and energy huq is emitted, an electron makes a 
transition from an initial state to a final state k' = k — q, whose rate is given 
by the Fermi golden rule: F^g^g, = ^\gq\'^6{Ej^ - hug - Ep)[nB + 1], where 
nB is the Bose distribution function for the phonons and gq is the coupling 
to be specified later. The rate of energy loss is 

= 2j2hujqT.,:^^,,[fiEk)il - f{E',))l (28) 

k,k' 

where f{Ek) is the Fermi distribution function for electrons, and the factor 
of 2 comes from the spin degeneracy. Following Ref. we replace the 

summations by appropriate integrals over k and k', and carry out first the 
integration perpendicular to the Fermi surface, yielding: 



lQT^^n\} J J [en^,/ksT. - 1] [1 - e-^^/'^BTq ' ^^^^ 



23 



where kf and vj are the electron wave number and speed at the Fermi surface, 
ds and ds' are elements of the Fermi surface. 



A similar consideration for the phonon absorption process yields an ex- 
pression for the rate of energy gain, which is the same as the above except 
that Tg and Ti have interchanged their places. The net rate of energy loss 
can then be written as Q = F{Te) — F{Ti), where 

We can replace one of the surface integrals by the surface integral of the 
relative variable k — k', then the other surface integral just yields a factor of 
the area of the Fermi surface. The relative surface integral can be further 
transformed to an integration over q, by using the relation qdq = kjsinOdO, 
where 6 is the angle between k and k' on the Fermi surface. The result is 

F(T ) - r .d. ^^^^^V rsii 



For simple estimates we may use Igg]"^ = |0 for the deformation 



potential coupling, where NZ is the total number of electrons in the metal. 
At low temperatures, we may further use Ug = ciq, and replace the upper 
limit of the integration by oo. The result is the famous law for the cooling 
rate: Q = J:V{T^ - ), where E = 24.9£5f|^, with n = NZ/V being the 
electron density. 

We are now in a position to make some simple estimates. For aluminum, 
we have Ej = 11.7 eV, n = 18.1 x 10^^/ m^, and q = 6420 m/s, yielding 
S = 0.1 nW/K^//im^. This compares fairly well with the measured value of 
S = 0.2 nW/K^/yum^ in Ref. [^, considering the crudeness of our model 
and the difficulty in the experiment. For copper, we have Ef = 7 eV, n = 
8.47 X 10^7 m^, and q = 5010 m/s, yielding S = 0.2 nW/K^/fim^. Although 
this predicts correctly a larger cooling rate for copper than for aluminum, the 
numerical value is about an order of magnitude smaller than the experimental 



result in Ref. [37 
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4.3 Cooling Rate of a Quantum Dot 

We now consider the cooling rate of a quantum dot made of a 2DEG in 
a semiconductor heterostructure. We assume that: (1) the dot is clean and 
small enough that impurity scattering may be neglected, and large enough 
that the discreteness of the energy levels may be ignored; (2) the electron gas 
form a degenerate Fermi sea, which is the case when the electron temperature 
Te is below the liquid helium temperature and the Fermi energy Ef is at or 
above a few meV. The cooling is achieved through coupling to phonons in 
the substrate, which is held at a lower temperature T;. 



Following the method of the last section, the rate of energy loss by emitting 
phonons is given by 

O V fdl' !dl ^^^<iyW .32^ 

where as before is the Fermi velocity, and gq is the coupling function. 
The lower dimensionality of the electron gas is reflected as follows: A is the 
area of the dot, dl and dV are line elements of the Fermi circle, and we have 
lost a pair of factors of However, the phonons are three dimensional, 
so gj| = k — k' , and we have the extra sum over q^. Again, we may replace 
the line integrals by 27rA;/ times / dq\\. Then, together with the sum, the 
latter integral may be written as ^ / qdq^ where is the thickness of the 
substrate in the z direction. Therefore 

Combining with a similar expression for the rate of energy gain due to phonon 
absorption, we can again write the net cooling rate as Q = FiT^) — F{Ti), 
where crA{T^ — Ti) at low temperatures, where the coefficient is given by 

- irf^ J ^^^e^.ABT_i- (34) 

For deformation potential couphng, we have ALz\gq\'^ — hq^D'^ / {2pijjq), 
yielding: 

AD'^kfiksTf /•°° x^dx 



^(^) - « 2/5 4 2 / ^-T' (35) 
Stt-^Ti pcfvj Jo — \ 



25 



at low temperatures. This gives Q = adAlT^—Ti"), with ad = 19 fW /K^/ /im^ 
for a 2DEG from GaAs heterostructure, where we have used the effective mass 
m = 0.067me and Fermi energy Ej = 7 meV, together with the material 
parameters for p, D, and q (see Section 3). 

For piezoelectric coupling, we have ALz\gq\'^ = he^cip^ /{2eoerq). This is to 
be screened by a factor of (gyas)^, which is |(gas)^ after directional average. 
This yields a similar form for the net cooling rate as Q = apA{T^ — T^'), 
where 

3C(5)fc/fc| g^eV 

with C being the Riemman Zeta function, and C(5) — 1.03693. This is 
evaluated as cTp = 10.2 fW/K^//im^, using the parameters mentioned in 
the last paragraph. The coefficient for the total cooling rate is a = aa + 
ap = 30 fW/K^//im^. This result agrees very well with the experimental 
result of 1^ for a high mobility sample at Tg > 0.12 K. However, at lower 
temperatures and for low mobility samples, the theoretical result is too low 
to explain the experiment. More theoretical work needs to be done to take 
the disorder effects into account. 



4.4 Scaling of Electron Temperature and Thermal Er- 
ror 

In this subsection, we will show how the electron temperature scales with 
the device operation frequency and size, and estimate its magnitude for dif- 
ferent types of devices under experimental conditions. We will assume that 
V < 7ee, where 7ee is the hot electron relaxation rate due to interaction with 
the background electrons, so that a meaningful electron temperature may 
be defined for the electron gas in a quantum dot or metal island, although 
it may be higher than the lattice temperature. The electron temperature is 
determined by the condition of power balance. The heating power due to hot 
electron (hole) deposition depends on the relative magnitude of u, 7ee, and 
7eph, where 'jeph is the relaxation rate of a tunneled- in electron via coupling 
to phonons. If •jeph < v < 7ee, phonon induced transitions can be neglected, 
and a hot electron relaxes by coupling to the feeding lead for a time l/7ee 
before being thermalized with the background electrons. \i u < jeph < 7ee, 
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the phonon channel is also involved in the relaxation before thermalization. 
If u < < ■Jeph, the extra energy of the tunneled-in electron is dissipated 
before heating takes place. For the sake of rough estimates, we will assume 
in the rest of this subsection that the energy relaxation of a hot electron due 
to tunneling and phonon emission may be neglected. The cooling power due 
to net phonon emission from the thermalized and hot electron gas has been 
studied in the previous subsections. 



Consider an electron turnstile made of metal islands with operating fre- 
quency z/, volume V, and Coulomb charging energy U. The heating power 
due to the deposition of hot electrons and holes in the island is ~ uU. Assum- 
ing that the lattice temperature T/ satisfies Tj^ <C T^, we obtain the cooling 
power as ~ HVT^, where S is a dissipation constant ~ 0.2 nWK~^/im" 



e 5 

[PH . In the steady state, we have Tg = (f^)^- Using the facts that U = 



1^)5. Using the facts that U = ^ 

1 14 

and C ~ \^3, we obtain that Te ~ z/sf/s. Hence, the thermal exponent that 
determines the thermal error is found as 

viTT.)'. (37) 



where r] = ( ^^^^^^ ) ^ is a dimensionless number, independent of the operating 
frequency and the volume of the island because U is inversely proportional 
to the linear size of the island. Using the known experimental values V = 
0.1 /im^ and C = 1 fF and 1/ = 10 MHz, we find U = 0.07 meV, r/ = 2, 
and Te = 80 mK. The electron temperature is much higher than the ambient 
temperature T = 10 mK. 



A similar result can be obtained for a quantum dot device, but the param- 
eter rj is given hj T] = Note that 7] now weakly depends on the 
linear size of the dot but not on the operating frequency. For an estimate, we 
use the data in Ref. §, where A = 0.5 /im^, C = 0.24 fF, yielding = 340 
mK at z/ = 10 MHz, which is much higher than the reported 10 mK ambient 
temperature, rj is found as ~ 1.9 for the given data. 



In an electron pump, the amount of heating in a cycle is not simply given 
by the charging energy because a tunneling stage is controlled by ramping a 
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gate voltage. For a 3-junction pump (See Figure 3), the heating energy can 
be derived as 

AE =^{64RtCu)^, (38) 
according to Eqs.(35) and (38) of Ref. [T^. Again, we write the thermal 



exponent as ^^^^ 



which weakly depends on u and Rt- For typical values of Rt = 100 kfl, 
C = 1 fF, and z/ = 10 MHz, we have r] ~ 1.3. 



In view of the above scaling results, we write the thermal error of these 
devices exp( — ^^^) in the more convenient form of lO~^*(^/'^^^^'\ where U 
is measured in units of meV, and u in lOMHz. The conversion relation is 
rj* = 3.261]. The thermal errors are then on the order of 10~^'^, 10~^'^, and 
]^g-6.75 j^Yie above typical examples respectively. These estimates seem to 
compare very well with the errors found in experiments, but the reader is 
warned not to take the numerical values too seriously because of the crude 
nature of our estimates. The scaling forms should be more generally valid. 



5 Conclusion 

In summary, we have presented a detailed study and overview of the 
nonequilibrium effects and self heating in single electron transfer devices 
which are primarily based on the Coulomb blockade effect. The Coulomb 
charging energy is taken to be much larger than the temperature (electronic 
or lattice), but the level spacing in the excitation spectrum is not necessarily 
so. We have considered three types of devices: quantum dot electron turn- 
stile, metal island electron turnstile, and metal island electron pump. In the 
turnstiles, a hot electron may be deposited in a quantum dot or metal island 
with an extra energy on the order of the Coulomb charging energy. In the 
electron pump, the extra energy of the hot electron depends on the rate of 
pumping. 
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We have considered three major channels of hot electron relaxation: (1) 
via tunneling, (2) emitting phonons, and (3) exciting background electrons. 
Three distinct regimes are identified: non-equilibrium (7ee,7ep?t < self 
heating or partial equilibrium {■jeph, ^ < lee), and equilibrium {jee, ^ < leph), 
see Figure 4 for an illustration. For the quantum dot device used in |], 
7eph is estimated to be about 10^ /s, while ■jee is on the order of 10^ /s, 
so that it lies in the self heating regime for practical operating freqencies. 
However, jee may be reduced by depleting the background electrons in the 
dot, rendering the system to be in the other regimes. For a metal island, 
both 'jeph and jee are found to be on the order of 10^ — 10^ /s, so for z/ > 10^ 
Hz the system is in the nonequilibrium regime, while for slower frequency 
the system lies near the boundary between the equilibrium and self heating 
regimes. 

In the nonequilibrium regime, a hot electron relaxes only through tunnel- 
ing back and forth to the feeding lead (or island). A set of characteristic rates 
are found for the hot electron relaxation. Those rates corresponding to the 
trapping levels below the Fermi energy of the feeding lead are exponentially 
small at low temperatures. This is understood as due to the Pauli exclu- 
sion in the lead, which is more and more effective at lower temperatures in 
blocking the tunneling-out process. A distributed set of characteristic rates 
give rise to an overall power law frequency dependence in the trapping error, 
which is closely related to the electron transfer error. For energy independent 
couplings between the levels and the lead, the power law is linear up to a 
logarithmic factor. If the higher levels couple more strongly to the lead, a 
sublinear power law is obtained. In the continuum limit, the tunneling-out 
process is nonoperative, and therefore a trapped hot electron is unable to 
relax. 

We then considered the relaxation due to phonon emission for a hot elec- 
tron in a quantum dot, and studied its effect on the characteristic rates due 
to tunneling. Perhaps nonsurprisingly, those slow characteristic rates are 
found to be pushed up to and above the relaxation rate 'jeph due to phonon 
emission. When the device operation frequency u is slower than jeph, the 
hot electron has enough time to reach equilibrium with the phonons (the 
ambient), with a residual trapping error exponentially small as a function of 
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V. For completeness, we have also considered a continuum model for metal 
island and large quantum dots. 



There is then the self heating regime, in which the hot electron quickly 
reaches equilibrium with the background electrons, but its extra energy is 
turned into electronic heat, elevating the electron temperature Tg above the 
lattice temperature Ti. We have rederived the formula Q = T,V(T^ — Ti) 
for the energy dissipation rate (cooling power) to the phonons from a hot 
electron gas in a metal island of volume V , and obtained a simple expression 
for the coefficient S. The coefficient is estimated to be 0.1 nW/K^//im^ for 
aluminum and 0.2 nW/K^//im'^ for copper. These results are in reasonable 
agreement with experiments [^]. A similar expression, Q = aA{T^ — Ti), 



has also been obtained for a quantum dot of area A, where the coefficient is 
estimated as 0.03 pW/K^//im^ for coupling to three dimensional phonons, 



with the result in fairly good agreement with the experimental finding of 
for a high mobility sample. 



The elevated electron temperature is obtained by balancing the cooling 
power with the heating power, and is found to be ~ 100 mK for a metal 
island and ~ 300 mK for a quantum dot under typical experimental condi- 
tions, while the ambient temperature is much lower. Different scaling be- 
haviors have been obtained for the thermal exponent f//(fcsT'e) as function 
of operation frequency and size for different type of devices. The estimated 
thermal errors are consistent with experimental results. 



There are two ways to solve the nonequilibrium and self heating problem. 
One way is to increase the rate of heat dissipation from the dot or metal 
island. This can be done to some extent by device scaling and by engineering 
the shape of the device structure in such a way that the dot area or island 
volume is maximized for a given capacitance. The other way is to reduce 
heat production. The electron pump is better than the turnstile precisely 
because of less heat production, especially so at low operation frequencies 
[0]. One may also consider making a quantum charge pump (or electron 
load- lock), which utilizes discrete electron levels to trap electrons, and has 
the advantage of producing essentially no heat in the dot 0. 
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Finally, we emphasize that although our study were focused on single elec- 
tron transfer devices, our results may also apply to other types of Coulomb 
blockade devices, such as electrometers and single electron transistors 



42]. The role of device frequency v is now played by the electron current 



divided by the electron charge. 
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Appendix: Non-Markovian Effects 



The trapping errors displayed in Figure 6 are very small. One wonders if 
the semiclassical rate equations (1) are accurate enough for the prediction of 
these small quantities. One possible correction to the semiclassical rate equa- 
tions comes from the non-Markovian effect. This question was investigated 
in detail in Ref. |2^, and here we give a short summary: 



Starting from the non-Markovian master equation (7) in Ref. 



set of differential-integral equations for the level occupation probabilities 
were derived. These equations were then approximated by the semiclassi- 
cal rate equations plus leading order corrections. The trapping errors were 
then calculated with the improved rate equations, with the results (un- 
der conditions of figure 6): 6.27 x 10^^ 7.71 x 10-3,3.27 x 10-^,4.35 x 
10-^ 1.73 X 10-^ 9.79 X 10"^ for operating frequencies v = 1 GHz, 100 MHz, 
50 MHz, 10 MHz, 1.25 MHz, and 0.125 MHz respectively. These are in 
good agreement with the corresponding semiclassical results 5.88 x 10"^, 
7.12 X 10-^ 2.97 X 10-^ 3.85 x 10-^ 1.69 x 10-^ and 9.37 x 10"^ 
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Figure Captions 



Fig. 1. Electron trapping in a quantum dot electron turnstile, which 

consists of a quantum dot (middle), two leads (with Fermi levels Ej' an Ej^), 
and energy barriers separating them. The Coulomb charging energy in the 
dot is denoted as U, and the level spacing as AE. With the left barrier 
lowered and under the biases shown in the figure, it is energetically favorable 
for one and only one electron to tunnel into the dot from the left lead, but 
the final state can be any one of the several available states in the dot. 

Fig. 2. The circuit for a metal island electron turnstile device. The metal 
islands lie in between the tunneling junctions of capacitances C. Cg and Vg 
are the gate capacitance and votage, respectively. 

Fig. 3. The circuit for a metal island electron pump. Each island is now 
separately gated. Electron tunneling through each junction is controled by 
ramping the bias accross it. 

Fig. 4. Three different regimes in the parameter space of jgph and 7ee, the 
relaxation rates of a hot electron through emitting phonons and exciting the 
background electrons. In the equilibrium regime, a trapped hot electron may 
achieve equilibrium with the background electrons and the lattice phonons 
before the trapping cycle is finished. In the self heating regime, the hot elec- 
trons may quickly thermalize with the background electrons, but with the 
latter heated to a temperature above the ambient lattice temperature. In 
the non-equilibrium regime, a trapped hot electron can relax only through 
tunneling back and forth to the feeding leed with a set of distributed char- 
acteristic rates, the lower part of which can be extremely slow due to Pauli 
exclusion. 

Fig. 5. Trapping probability P{t) vs. ut (solid line). Also shown are 
the currents II (dot-dashed line) and Ir (dashed line) in units of u. The 
arrows indicate the times at which P{t) achieves its minimum, half value, 
and maximum. 
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Fig. 6. The solid curve is the approximate theoretical prediction from Eq. 
(3), which fits very well the results of an accurate numerical similation (di- 
amonds). The exponentials indicate the behaviors of the downward curving 
segments, with z/2 = 14 x 10^ /s, = 1.2 x 10^ /s, and 2/4 = 96 x 10^ /s. 
The overall behavior is power law like, with 1 — Max{P) ~ 
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Table Caption 

Table 1. Two sets of egienvalues Xj in the absence and presence of phonon 
induced transitions. The parameters taken are N — 10, 6 — 10, = 4, 
and rj^rf^ 10^0 /s. 
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